Three dimensionality in quasi-two dimensional flows: recirculations and Barrel effects 
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A scenario is put forward for the appearance of three-dimensionality both in quasi-2D rotating 
flows and quasi-2D magnetohydrodynamic (MHD) flows. We show that 3D recirculating flows and 
currents originate in wall boundary layers and that, unlike in ordinary hydrodynamic flows, they 
cannot be ignited by confinement alone. They also induce a second form of three-dimensionality 
with quadratic variations of velocities and current across the channel. This scenario explains both 
the common tendency of these flows to two-dimensionality and the mechanisms of the recirculations 
through a single formal analogy covering a wide class of flow including rotating and MHD flows. 
These trans-disciplinary effects are thus active in atmospheres, oceans or the cooling blankets of 
nuclear fusion reactors. 

PACS numbers: 47.32.-y, 47.65.-d, 47.27.nd 



I. INTRODUCTION 

Rapidly rotating flows and electrically conducting 
flows in homogeneous static magnetic fields share a 
remarkable property: both tend to two-dimensionality. 
In the former, inertial waves propagation promotes in- 
variance along the rotation axis In the latter, when 
magnetic field perturbations induced by flow motion are 
negligible (in the quasi-static MagnetoHydroDynamic 
(MHD) approximation), electric eddy currents damp 
velocity variations along the magnetic field lines (l4j . 
This feature is crucial because 2D and 3D flows have 
radically different dissipative and transport properties, 
especially when turbulent. While in 3D turbulence, 
energy follows a direct cascade from large to small 
structures where it is dissipated by viscous friction, 
2D turbulence proceeds through an inverse cascade 
that accumulates energy in large structures where it 
is dissipated by friction on the boundaries [H, [2(|. In 
planar fluid layers, this mechanism is suppressed as 
three-dimensionality emerges [l6[ and its suppression or 
appearance are respectively determined by the presence 
or absence of fluid motion along the third component 
[H, |22| . The a priori distinct questions of the number of 
components of velocity field (2C/3C) and the of spatial 
directions in which it varies are thus tightly linked and 
determine whether flows obey 2D or 3D dynamics, how 
and when they may switch between them. 
The fundamental difference between 2D and 3D states 
places these questions at the centre a vast array of 
problems across disciplines: atmospheric and oceanic 
flows, electromagnetic flow control in metallurgical 
processes, even the dynamo problem. Despite their 
importance, the transition mechanisms between 2D and 
3D states, in rotating, MHD or other flows are not 
understood yet, and an important gap exists between 
theory and experiments. Numerical simulations show 
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that with periodic or free-slip boundaries, instabilities 
can break down strictly 2D MHD structures and lead to 
2D-3D intermittency [j,|2l|. In rotating flows, non-linear 
transfer occurs from 3D to strictly 2D modes [l7j]. In 
experiments, shallow layers cannot be strictly 2D but 
only quasi-2D, either because of the three-dimensionality 
introduced by viscous boundary layers or because of 
variations of the same order in the core flow: firstly, such 
quasi-2D flows ignite 3D recirculations, that are crucial 
to the 2D-3D transition because of their influence on the 
inverse cascade and also because they nurture small-scale 
3D turbulence [ljj. Secondly, in MHD channel flows, 
eddy currents between the boundary layer and the core 
were shown to induce a quadratic three-dimensionality 
lol Il2| , an effect recently observed experimentally, too 
8j. Currently, the mechanisms that promote one or the 
other of these effects are unknown. Viscous boundary 
layers most certainly play a role in promoting 3D 
recirculations in decaying and steady vortical flows [Tfjj . 
with and without background rotation Nevertheless, 
in the absence of Coriolis or Lorentz force, confinement 
alone suffices to trigger them Q, even in the absence of 
boundary layer friction. These results stress the need for 
a clarification of the roles both of external forces such as 
the Lorentz or the Coriolis force on the one hand, and of 
confining boundaries on the other hand. 
In this Letter, we examine the simple configuration of a 
symmetric, plane channel in a transverse field, bounded 
by two no-slip walls distant of 2H (fig. [T]), to reveal 
how walls induce three-dimensionality in quasi-2D flows. 
The governing equations are written in a general form 
that defines a large class of flow encompassing not only 
quasi-static MHD and rotating flows, but also rotating- 
MHD flows relevant to geophysical dynamos. In doing 
so, a formal analogy is defined within this class of flows, 
which share a common tendency to two-dimensionality. 
A hierarchy of mechanisms is singled out that both 
explains how 3-Component (3C) motions are ignited 
and how three-dimensionality appears as a variation of 
physical quantities across the channel. 
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2. 3C velocity >~ 

3. 3D current y / 




FIG. 1. Channel flow configuration. The three steps of the 
first Barrel effect are schematically represented. In the second 
type of Barrel effect, active in MHD flows, the steps are iden- 
tical with the roles of the current and the velocity swapped. 



II. GOVERNING EQUATIONS 

Following the formalism of rotating and quasi-static 
MHD flows, we first define a static homogeneous field 
F = Fe z (vertical for convenience), orthogonal to the 
channel, that pervades the fluid domain and interacts 
with a current c to exert a force c x F on the flow. In 
rotating flows (resp. quasi-static MHD), F is twice the 
background rotation vector (resp. the magnetic field) 
and c is the flow current (resp. electric current), that 
appear in the Coriolis (resp. Lorentz) force. The current 
is linked to the velocity field u by a phenomenological law 
such as Ohm's law in MHD, that takes the general form 
c = J-(u) . The corresponding non-dimensional equations 
for an incompressible fluid (density p, viscosity v) depend 
on two dimensionless groups, built on H, the reference 
velocity U and reference current C: e = pU 2 / (HCF) and 
e 2 = pvU I (H 2 CF) respectively express the ratio of iner- 
tia and viscous forces to c x F, and both are smaller than 
unity. This class of analogy is summarised in table HI Wc 
focus on regimes where deviations to the leading order 
play a key role, and in particular those due to inertia, 
i.e typically e < 1 and 10~ 2 < t v << 1. Denoting any 
quantity g in the core (i.e. outside the wall boundary 
layers) as <?, the equations write: 

e—u x + V_u5 = c ± x e z + t 2 v (d 2 zz + V 2 ± )u x , (1) 



d z ii z = -V ± • Uj_, 
d z c z = -V ± • c ± , 
c = F(u). 



(2) 

(3) 
(4) 
(5) 



Vector quantities and operators are split into their com- 
ponents along e z (subscript z) and in the (x, y) plane 



(subscript _L). In rotating flows, c disappears from the 
governing equations, as ([5]) simply becomes c = u. In 
low-i?m MHD, on the other hand, ((5]) represents Ohm's 
law: 



(6) 



where the core electric potential <f> is usually obtained as 
solution of the Poisson equation that follows from sub- 
stituting © into Q P3. 



III. TENDENCY TO TWO-DIMENSIONALITY 
IN THE CORE FLOW 

We shall first clarify the interplay between the respec- 
tive dynamics of (c_l,c z ) and (u±,u z ) in the core. At 
the leading order (e°, e°) when e -> 0, e v -> 0, dJ) and (J2J 
readily imply that Cj^ and p are independent of z (i.e. 
2D). Denoting leading order quantities with the super- 
script g°, 



c () L =e z x \7 ± p°, 



d z p° = 0. 



(7) 
(8) 



The problem symmetry and (j4]) further imply that the 
current is exclusively horizontal: 



0. 



(9) 



Two-dimensionality of the core velocity isn't ensured at 
this point, but depends on the nature of §5§. For the 
time being, we shall assume that z-dependence doesn't 
explicitly appear in T . In this case, is indeed 2D. 
When applied to rotating and quasi-static MHD flows, 
this result recovers the property that both are quasi-2D 
in the absence of inertial and viscous effects in the core. 
The z— dependence of c± appears from (p} and ([2|): 



d z 



dt 



elV 2 



(10) 



Since the leading order core velocity is 2D, cj^ must in 
fact be 2D at least up to 0(e,e 2 ), so that using (0} and 
the problem symmetry yields: 



-zc z (-l) + 0(e,el). 



(11) 



Similarly, the problem symmetry and ([3]) imply that 



-zu z (-l) + 0(e,el). 



(12) 



From (|10p and (fTTj) . the appearance of z— variations 
in cj^ and c z is determined by the flow and the current 
injected into the core from the boundary layers that 
develop along the walls u z (— 1) and c z (—l). Only then, 
can z— variations in appear through (|5|), at the same 
order as c^. Therefore, we now need to analyse the wall 
boundary layers. 
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Field F 


current c = J-(u) 


e 




Wall boundary layer (u l j_ =) 


Quasi- 
static 
MHD 


magnetic 
field 
B 


electric current 
j = a{-V(f) + u x B) 


Stuart number 
N ~ e" 1 - ° B * H 

pU 


Hartmann 

number 
Ha = e- 1 = 

V pv 


Hartmann layer [9j 
u° x (-l)(l-e-f) 


Rotating 
flows 


double 
rotation 

2to 


flow current 
pu 


Rossby number 
fib - f - -V— 

™ — t — 2HS2 


Ekman number 
E - <? - K 


Ekman layer [rf]| 

ui(-l)(l-e"^cos^) + 
e z x Uj_(— l)e~75 sin 



TABLE I. Analogy table giving the expressions of generic dimensional quantities F, c, and non-dimensional numbers e„ and e 
for quasi-static MHD and rotating flows, a and <j> are the electric conductivity of the fluid and the electric potential. 



IV. RECIRCULATIONS DRIVEN IN THE 
BOUNDARY LAYER 

In the vicinity of walls, horizontal viscous friction must 
balance c x F to achieve the no slip boundary condition, 
and this imposes the thickness of the resulting boundary 
layer to scale as e v (respectively the Ekman and Hart- 
mann layers in rotating and MHD flows). The govern- 
ing equations are accordingly rewritten in the bound- 
ary layer near the wall z = — 1 using stretched variable 
( = e,7 1 (z + 1), to reflect that viscous friction becomes 
0(1). Denoting any quantity g in this region as g, 

+ V_lj3 = c± x e z + (d 2 cc + e 2 v V 2 x )u± , (13) 



ej t u z + d z p=(d 2 (c +elVl)u z , (14) 

d(U z = -£„V_l • u_l, (15) 

9qc z = -e„V_L • c_l, (16) 

c = ^(u). (17) 

The no slip condition at the wall is written 

u x (-l) = 0, (18) 



while the boundary conditions for c and u z shall be left 
unspecified for the time being, thus allowing for possi- 
ble wall-injection of current or mass into the fluid. Core 
and boundary layer variables must also satisfy a match- 
ing condition. Since they are not explicitly expanded in 
powers of e, e„, the general matching condition for any 
quantity g put forward by 0] simplifies to [l3j : 

Urn J(z+1)= Urn g(0- (19) 

c z (— 1) is obtained by integrating \7± x (Tl~3|) using (fTS)) . 
(fT9| . noting that d^u> z (— 1) = e u d z Cu z {— 1): 

c z {-\) = cf + e„9 cW f + 0(66,, e 2 v ), (20) 

where the superscript W refers to values taken at the 
wall. The current injected at the wall would be 
determined by a boundary condition for the current, 
here unspecified. Eqs. (|2"0"|) and (fTTj) reveal the two 



main mechanisms that can feed eddy currents in the 
core: injection of current at the wall and rotational wall 
friction contribute respectively at the leading order, and 
at 0(e„). Importantly, (jTOJ) implies that although these 
recirculations make c 3C, they do not directly affect 
the z— dependence of c^(z), which remains 2D up to 
0(e, el) at least. 

Strikingly, if cf ^ 0, (20]) contradicts ©: a linearly 
z-dependent vertical current in the core leads to di- 
verging horizontal currents that induce a rotational 
force c± x e z , which cannot be balanced by the sole 
leading order pressure gradient in (|7j). To resolve this 
paradox, viscous or inertial effects must exist in the core 
to oppose c± x e 2 . These can therefore not be neglected 
and, in fact, prevent two-dimensionality at the leading 
order. Quasi-static MHD provides a well understood 
manifestation of this effect: MHD flows are often driven 
by injecting electric current through point-electrodes 
embedded in an otherwise electrically conducting wall 
[l8| . Above such an electrode develops a vortex of 
rotation axis e z , with a viscous core that is 3D at the 
leading order. Injecting current at the wall therefore 
directly prevents quasi-two dimensionality, at least 
locally. Having now singled out this important effect, we 
shall assume = 0, unless otherwise specified, for the 
reminder of this Letter and focus on higher order effects. 
At order 0(e v ). (f2"0")) generalises two classic properties of 
Ekman and Hartmann boundary layers: in rotating flows, 
friction in the Ekman layer gives rise to secondary flows 
in the core u z {— 1) ~ ^E 1 / 2 ^^— 1), by Ekman pumping 
[lH . In quasi-static MHD, c\u^ = u±(-l) + (D(e v , e), so 
([2"D")) expresses that vorticity in the core drives an electric 
current 7 2 (—l) — i?a _1 w 2 (— 1) out of the Hartmann 
layers [9| (these two results are recovered using the 
leading order solutions of (JT5JT7J), u±, given in table |T]). 

We shall now turn our attention to the determination 
of u z (— 1), which controls d z c±(z). First, from (fT4| and 
(fT"5)) . the pressure is constant across the boundary layer: 

V 1 _p{Q = V ± P{-l) + 0{ee v ,el). (21) 

The horizontal pressure gradient in the boundary layer 
thus results from the balance of forces in the core. (f2"Tj) 
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and (fT3")) express that each of these forces alters the local 
balance between viscous forces and cxe 2 in the bound- 
ary layer, and thereby drives horizontal jets. Should 
these jets be horizontally divergent, they in turn induce 
a vertical flow from the boundary layer to the core. The 
equation for u z thus follows from V i - (|13|l and (fT5|) . and 
u z (— 1) is obtained by integration, using (TT^J) and ([12")) : 

(1- e v =uY + e u u° + e u eul + 0(4), (22) 

p OO p OO P oo 

= / / / Vix[ci- M-l)] ■ e,dCd&d&, 
Jo Jc Jc 2 



C •'G 
OO />00 /»oo 



^C: 



V, • 



d - ( -n d - 

dt^-v - dt ul - 



dCdCidd 



Eq. (|22p singles out three possible origins of secondary 
flows between core and boundary layer: the flow di- 
rectly injected at the wall doesn't interact with 
the boundary layer and is integrally transmitted to 
the core where it affects the flow at the leading order. 
The symmetric boundary conditions at the walls imply 
through (fl"2"|) that it creates two symmetric recirculations 
there. Should no flow be injected at the wall (u™ = 0), 
secondary flows can result from a rotational current 
as u z {— 1) = t v v!t + 0(ee y , el). Unlike wall-injected 
flow, the corresponding vertical flow builds up from 
horizontally divergent jets in the boundary layer but 
recirculates in the core in the same way. It leads to 
secondary flows at order e v there, a scaling they inherit 
from the thickness of the boundary layer where they are 
created. Finally, if the core current is curl-free at the 
leading order, then = 0. Inertial effects of order e 
then take over as the dominant mechanism that drives 
jets in the boundary layers, resulting in a secondary flow 
iiz(— 1) = ee v u z + 0(ef,). Importantly, no recirculation 
occurs if uY — Uz = u{ = 0, so in contrast to quasi-2D 
flows not subject to an external homogeneous field F Q , 
they are not triggered by confinement alone. 
In rotating flows, c = u and V x = u> z e z ^ 
in general (see table fl}: the mechanism responsible 
for secondary flows is the Ekman pumping men- 
tioned earlier when analysing eddy currents fed by 
c z (— 1). Still, application of (|2"2"|) again recovers the 
well-known result that u z {— 1) = e v u^ + 0{ee Vl e 2 v ) — 

^E x l 2 u z {-\) + (ibEVa^pl. 

In quasi-static MHD, by contrast, = so in the 
inertialess theory of the Hartmann layer, no flow escapes 
to the core @. When inertia is taken into account 
though, ([22j) yields u z (-l) = ee v u\ + Q(ee Vl e 3 „) = 
-(5/6)Ha~ 1 N- 1 V • [<(-!) ■ Vu° (-1)" 



V. THE "BARREL" EFFECTS 

Having expressed the current and mass flow that feed 
into the core, we are now in position to return to the 
core flow equations analysed at the beginning of this 



Letter and determine the conditions of appearance of 
z— dependence in core quantities. First, since c_|_ appears 
at a higher order than uj^ in (JTJI3J) , and since we further 
assumed that z— dependence didn't appear explicitly in 
JF, z-dependence cannot appear at lower order in than 
in cj^. Consequently, at the first order at which it ap- 
pears, the second term in the expression of d z c± (fTU)) 
vanishes. By virtue of the symmetric boundary condi- 
tions, and depending on (|22[) . three-dimensionality may 
appear under either of the three forms: 



(23) 



c ± = fi ± (0) + e l-|vx«r + 0(ee v , e 2 u ), 
c x = c x (0) + ce^lvxit? + 0(e 2 e u ,el), (24) 

c ± = MO) + eh^jVM + O(el). (25) 

Physically, these equations express that the z-linear 
vertical flow created in the core by the mass flow 
ejected from the boundary layers builds up a z-quadratic 
pressure in the core. A quadratic current must in turn 
be drawn in the core, for the force cx X e z to be able 
to balance the corresponding quadratic component of 
the horizontal pressure gradient. The z-dependence of 
is determined by the nature of J-: in rotating flows 
for example, = .F(u_l) = u±. Consequently, uj_ 
also depends quadratically on z, because does. Core 
structures are thus not columnar as usually assumed in 
quasi-2D flows but rather 6arre/-shaped. 
The Barrel effect has not yet been observed as such in 
rotating flows, but was discovered in the MHD numerical 
simulations of [Io| and predicted theoretically by (l2l |. 
We shall now show that the same general mechanism 
is at play in both cases. To this end, we must invoke 
a second analogy between MHD and rotating flows 
and release the assumption that J- doesn't explicitly 
depend on z: in MHD, the non-dimensional expression 
of T is given by Ohm's law c = j = — V0 + u x e z 
(table [I]). At the leading order, z-dependence is still 
absent in J- . so all core quantities remain 2D. Appli- 
cation of (TlTj) and ([2"U)) provide the expression of the 
vertical current there: c z = _ffa _1 zw° + 0(Ha~~ 2 , A -1 ). 
According to Ohm's law, this vertical electric current 
induces a quadratic component of the electric potential 
4> = {z 2 /2)Ha- 1 Lb° z + 0(Ha~ 2 , N^ 1 ), which introduces 
an explicit z-dependence in T at this order. Since, 
however, (fl"U|) still implies that must be 2D at 
C(i/a _1 ), Ohm's law demands that be quadratic 
at C(i/a _1 ). This mechanism is analogous to the 
Barrel effect identified previously, with j z , (f> and Ohm's 
law taking over the respective roles of u z , p and (|2|). 
The first Barrel effect was driven by vertical flows out 
of the boundary layer into the core, induced by 2D 
rotational currents. The MHD Barrel effect, by contrast, 
is driven by vertical currents out of the boundary layer, 
due to 2D vorticity. This second analogy underlines 
that although one common formalism explains why 
rotating and MHD flows are either 2D, quasi-2D or 
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3C, the origins of their "Barrel" three-dimensionality 
are still formally analogous but involve a distinct analogy. 



VI. INFLUENCE OF THE BOUNDARY 
CONDITIONS 

Up to this point, the analysis has been confined to 
a symmetric channel flow bounded by no-slip walls. In 
the quest for a flow closer to strict two-dimensionality, 
however, these conditions were often altered in a number 
of experimental and numerical studies of quasi-2D flows. 
We shall now examine how the driving mechanism for 
3D recirculations and barrel effects are affected by these 
different boundary conditions in two main cases: 



A. No slip bottom wall and upper free surface 

The upper wall at z = 1 is replaced with a non- 
deformable free surface (i.e. it remains flat). This config- 
uration was studied experimentally and numerically by 
The author showed in particular that free surface 
deformation didn't play an important role in the appear- 
ance of three-dimensionality, thus justifying the relevance 
of this boundary condition. Further assuming that no 
current escapes through the free surface, the boundary 
conditions at z = 1 become: 

a a uj_(l)=0 u»(l) = c„(l)=0. (26) 

In the channel with upper and lower no-slip walls, the 
problem symmetry implied d z u±(Q) = 0, u z (0) = 0, and 
c z (0) = 0, so the flow with a non-deformable free surface 
is simply that in the lower half of the channel found pre- 
viously and the physical mechanisms remain the same. 



B. Upper and a lower non-deformable free surfaces 

This second, more ideal variant was investigated nu- 
merically with and without background rotation [l|, Q . 
(jT5J) is replaced with 

0,ul(-1) = O. (27) 

The direct consequence is that the viscous boundary layer 
at z = —1 is no longer present, so this boundary condi- 
tion applies directly to the core flow (i.e. u = u and 
c = c). Again, the boundary condition at z = —1 for 
u z and c z can be left unspecified, to allow for possi- 
ble injection of mass and current through the bottom 
free surface (The injection can be made symmetric as for 
the channel flow, by relaxing the condition c z (l) = lo- 
cally) . This configuration is also symmetric with respect 
to the z = plane so equations (fTTj) and (fT2"]) for the core 
remain valid. From pop , injecting current at z = — 1 



therefore still leads to current recirculations and to three- 
dimensionality, both at the leading order, and locally pre- 
vents quasi-two dimensionality. Injecting mass also feeds 
recirculations at the leading order, but incurs a barrel 
effect at 0(e), also from (TO]). If u z (-l) = c z (-l) = 0, 
it follows from (TlTj) and (|12p that c z = u z = 0: both 
vertical flows and vertical current become severed from 
their only source in (fTTj) and ([T2"]) , so no recirculation can 
occur, and the flow remains 2C. Since the barrel effects 
are driven by these recirculations, they are shut down 
too and the flow is indeed strictly 2D. Such conditions 
are ideal but they reveal the crucial role played by the 
boundary layers as sources of three-dimensionality in the 
presence of an external field F. This result comes in stark 
contrast to similar configurations without external field 
where confinement alone incurs three-dimensionality Q. 
In practice, Q proposed to reduce the influence of the 
boundary layer by resting the active fluid layer on an- 
other one, acting as buffer. The interface between the 
two fluids was assumed non-deformable but the bound- 
ary condition for the bottom boundary of the active layer 
wasn't exactly that of a free surface. Local viscous fric- 
tion still existed and a viscous boundary layer developed 
there, albeit incurring less lower friction than if a solid 
wall were present. The authors noticed reduced but per- 
sistent three-dimensionality. Our analysis shows that 
in the presence of an external field, the same effect is 
present, as reducing the friction at the bottom boundary 
in (|20p would immediately damp both the current recir- 
culations through (|TTj) and the subsequent barrel effect. 



VII. CONCLUSION 

In this Letter, a quantitative analysis of the condi- 
tions of appearance of three-dimensionality in quasi-2D 
flows under the influence of an external field was con- 
ducted. It was shown that in the absence of any three- 
dimensional forces, three-dimensionality appears under 
the two forms of 3D recirculations and transverse vari- 
ation of fluid quantities. We called the latter "Barrel 
effect" and proved that it was driven by the former. Al- 
though a second order correction in the asymptotic sense, 
the Barrel effect becomes crucial at the transition be- 
tween 2D and 3D regimes. In turbulent flows, three- 
dimensionality is indeed driven by inertia and therefore 
appears when e becomes closer to unity. This implies 
that three-dimensional instabilities that lead to a fully 
3D state may not develop on a base flow made of colum- 
nar structures, as in strictly 2D flows, but rather on a 
flow with a quadratic profile that cannot be neglected 
and whose stability properties are currently unknown. 
Finally, this Letter stresses that in quasi-2D fluid lay- 
ers pervaded by a uniform transverse field F but sub- 
ject to no explicitly 3D force, 3D recirculations and flow 
three-dimensionality are linked and occur exclusively be- 
cause of the presence of boundary layers along the con- 
fining walls. This complements the findings of [l|,Q, who 
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showed that in the same configuration, with and with- 
out background rotation, three-dimensionality could still 
be observed in the absence of boundary layer friction ei- 
ther if the forcing was itself 3D, or without background 
rotation, because of confinement alone. 
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